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Abstract. We establish the time decay rates of the solution to the Cauchy problem for the 
compressible Navier-Stokes-Poisson system via a refined pure energy method. In particular, 
the optimal decay rates of the higher-order spatial derivatives of the solution are obtained. The 
H~''{0 < s < 3/2) negative Sobolev norms are shown to be preserved along time evolution 
and enhance the decay rates. As a corollary, we also obtain the usual L^-L^{1 < p < 2) 
type of the optimal decay rates. Compared to the compressible Navier-Stokes system and 
the compressible irrotational Euler-Poisson system, our results imply that both the dispersion 
effect of the electric field and the viscous dissipation contribute to enhance the decay rate of 
, the density. Our proof is based on a family of scaled energy estimates with minimum derivative 

CNj ■ counts and interpolations among them without linear decay analysis. 
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1. Introduction 



The dynamic of charged particles of one carrier type (e.g., electrons) in the absence of mag- 
netic effects can be described by the compressible (unipolar) Navier-Stokes-Poisson equations 
Oh ; (NSP): 

dtp + div(pu) = 

dt{pu) + div(/?n ®u) + Vp{p) — pAu — {p + A)Vdivu = pV$ , . 

A<^ = p-p 
^ {p,u)\t=o = {po.uo). 

Here pit, x) > 0, u{t, x) represent the density and the velocity functions of the electrons respec- 
tively, at time t > and position x G M'^. The self-consistent electric potential ^ = ^{t,x) 
^ ■ is coupled with the density through the Poisson equation. The pressure p = p{p) is a smooth 

, function with p'{p) > for /? > 0. We assume that the constant viscosity coefficients p and A 

^ I satisfy the usual physical conditions 

P>0, X+Ip>0. (1.2) 

In the motion of the fluid, due to the greater inertia the ions merely provide a constant charged 
background p > 0. For simplicity, we take p = 1 and assume that p'{l) = 1. 

The main purpose of this paper is to investigate the influence of the electric field on the time 
decay rates of the solution compared to the compressible Navier-Stokes equations (NS). We first 
^ ■ review some previous works on the global existence of the solutions to the NSP system. For the 

pressure p{p) = p^ with 7 > 3/2, the global existence of weak solutions was obtained by [19] 
when the spatial dimension is three. Later this result was extended by [31] to the case 7 > 1 
when the dimension is two, where the authors introduced an idea to overcome the new difficulty 
caused by that the Poisson term pS/^ may not be integrable when the 7 is close to one. The 
constraint of 7 is somewhat optimal in the sense of the well-known framework of weak solutions 
to the NS system [UJI^n]. The global existence of small strong solutions in Sobolev spaces 
was shown in [22] in the framework of Matsumura and Nishida [25j, while global existence of 
small solutions in some Besov spaces was obtained in |1UI W2\ . 

The convergence rate of the solutions towards the steady state has been an important prob- 
lem in the PDF theory. The decay rate of solutions to the NS system has been investigated 
extensively since the works [251 [26l III] , see for instance [MllMlliailHllITlliaillllllSllIlllllIIl 
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\T2\ [21} I34| and the references therein. When the initial perturbation pQ — 1,uq G L''^ n with 
p € [1,2] (Indeed, in those references p is near 1 and A'^ > 3 is a large enough integer for the 
nonhnear system.), the optimal decay rate of the solution to the NS system is 



_3 1 



\\{p-l,n){t)\\^,<il + t)-n-r-V. (1.3) 

Recently, the decay rate of solutions to the NSP system was investigated in [22l |36l [331 IE] • It is 
observed that the electric field has significant effects on the large time behavior of the solution. 
When the initial perturbation po — l,uo G D with p G [1,2], then the optimal decay 
rate of the solution to the NSP system is 

\\ip-m)\\L^<a + tr'^^~''^'^ and ||n(t)||^2<(l + t)"^(^^)+i (1.4) 

This implies that the presence of the electric field slows down the decay rate of the velocity of 
the NSP system with the factor 1/2 compared to the NS system. The proof is based on that 
the NSP system can be transformed into the NS system with a non-local force term 

dtp + div(/?n) = 

dtipu) + div{pu (^u) + Vp{p) - pAu - {p + A)Vdivu = pVA'^ip - 1) (1.5) 
ip,u)\t=o = ipo,uo). 

By the detailed analysis of the Fourier transform of the Green function for the linear homoge- 
neous system of (jl.Sp . we may have the following approximation for the Fourier transform of 
the solution, by refining the estimates (3.3)-(3.4) of [22j, 

/ 0(l)e~(^+^^)l«l^* (l^ol + iCllnol) , lel < V, 



O(l)e-^°*(|^o| + !^o|), lel >^?, 



and 



/ o(i)e-'^i«i*(ier^i^oi + |noi), iei<r?, 

u(t,t) ~ < „ (1.7) 

1 O(l)e-^"*(|^o| + |^o|), \i\>V- 
Hereafter we may sometimes write g = p — 1. i?o>Oisa constant and r/ > is a small but 
fixed constant. Then the linear optimal decay rate (|1.4p follows if go,UQ G with p € [1, 2]. 

However, in this paper we will give a different (contrary) comprehension of the effect of the 
electric field on the time decay rates of the solution. The key motivation is that if we take p = 2 
in the time decay rate ()1.4p , then we should get that the norm of u grows in time at the rate 
(1 + t)^/^! This seems unsuitable since the NSP system is a dissipative system. The reason why 
this happened is that to derive (|1.4p with p = 2 it only assume that Qo,uo € L^, but from the 
point of view of the energy structure of the NSP system it is natural to assume that V$o ^ L,^- 
The linear energy identity of the perturbation form of (jl.ip reads as 

~[ |p- Ip + [lip + |V^'pdx+ / p\Vu\'^ + {p + X)\divufdx = 0. (1.8) 
2 at J^3 JR3 

By the Poisson equation, the condition V^>o € is equivalent to that A^^^o S L^- Motivated 
by this, we instead assume that A~^^>o, ""0 ^ with p G [1, 2], then by ()1.6p ~ ()1.7p . we have the 
following optimal decay rates for the linear NSP system: 



3/ll\l 3/11 



+ ^ and \\u{t)\\^, < (1 + 1) Hp V. (1.9) 

In this sense, the electric field does not slow down but rather enhances the time decay rate of the 
density with the factor 1/2! This can be understood well from the physical point of view since 
we get an additional dispersive effect from the repulsive electric force. This is also consistent 
with [7] in the study of the compressible Euler-Poisson equations. 

In the usual L^-L^ approach of studying the optimal decay rates of the solutions, it is difficult 
to show that the norm of the solution can be preserved along time evolution. Motivated 
by [8], using a negative Sobolev space {s > 0) to replace L'p space, we developed in [9] a 
general energy method of using a family of scaled energy estimates with minimum derivative 
counts and interpolations among them (but without linear decay analysis) to prove the optimal 
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decay rate of the dissipative equations in the whole space. An important feature is that the 
norm of the solution is preserved along time evolution. The method was applied to classical 
examples in [9] such as the heat equation, the compressible Navier-Stokes equations and the 
Boltzmann equation. In this paper, we will apply this energy method to prove the optimal 
decay rate of the solution to the NSP system p.ip . 

Notation. In this paper, with an integer £ > stands for the usual any spatial derivatives 
of order £. When ^ < or £ is not a positive integer, stands for defined by (|A.9p . 
We use H^{M.^),s € M to denote the homogeneous Sobolev spaces on M'^ with norm 
defined by (jA.lOp . and we use H'^{M.^) to denote the usual Sobolev spaces with norm \\-\\jjs and 
L^'(M^), 1 < p < oo to denote the usual spaces with norm IHI^^p- We will employ the notation 
a < 6 to mean that a < Cb for a universal constant C > that only depends on the parameters 
coming from the problem, and the indexes and s coming from the regularity on the data. 
We also use Co for a positive constant depending additionally on the initial data. 
Our main results are stated in the following theorem. 

Theorem 1.1. Assume that po — 1, uq, V<^o ^ for an integer N > 3 and 

(1.10) 



{po — I) dx = (neutrality). 

i 

Then there exists a constant 6o such that if 

Wpo - M\h» + \\uo\\h3 + llV^olli^s < So, (1.11) 
then the problem p.ip admits a unique global solution (p, u, V$) satisfying that for all t > 0, 



||(p-l)(t)||^. + ||n 
< C 



+ ||V$ 



UN + 



f \\{P - l)(r)||^iv + ||Vn(r)||^^ + ||VVcI>(r)||^^ dr 
Jo 



\Po 



If further, pQ — 1, uq, V^q H ^ for some s € [0, 3/2), then for all t >0, 



wip-mnl^ 

and the following decay results hold: 
+ V^n(t) 



+ u 



+ ||V$ 



H 



v'{p-m) 

and 



+ 



VV$(t) 



<Co(l + t)~^ 



for i = 0,. 



V\p-l)it) 



L2 



<co(i+r^ 



for £ = 0, 



A 



(1.12) 
(1.13) 

,A-1, 
(1.14) 

(1.15) 



Note that the Hardy-Littlewood-Sobolev theorem (cf. Lemma [A. 5 P implies that for p E (1, 2], 
L'P C H~''^ with s = 3(^ ~ ^) ^ [0' 3/2). Then by Theorem I l.H we have the following corollary 
of the usual U'-L? type of the optimal decay results: 

Corollary 1.2. Under the assumptions of Theorem \1.1\ except that we replace the as- 
sumption by that po ~ ^iUqjV^o € for some p € (1,2], then the following decay results 
hold: 

V'{p-l){t) 



+ 



V'u{t) 



+ 



W<^{t) 



HN+l- 



<Co{l+t) 



and 



V\p-l){t) < Co{l + t)-^''p-'+^2) for£ = 0,...,N 



for. 



2. 



0, 



Here the number ap^i is defined by 



3 /I 



A-1, 
(1.16) 

(1.17) 



(1.18) 
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The followings are several remarks for Theorem 11.11 and Corollary ll.2[ 
Remark 1.3. We remark that by the Poisson equation, 



1 



+ ||A Vo 



S > 0, 



and 



\Po 



11 



+ ||A Vo 



I LP 



[1,2]. 



(1.19) 
(1.20) 



Compared to the study of the NS system, such norms on A^^qq are additionally required. But 
they can be achieved by the natural neutral condition (|1.10p . see the proof in pp. 263-264 of [7]. 

Remark 1.4. Notice that for the global existence of the solution we only assume that the 
norm of initial data is small, while the higher-order Sobolev norms can be arbitrarily large. This 
is an improvement of [9] where we required the smallness of the H^^^^"^ norm of initial data. 
Also notice that we do not assume that or norm of initial data is small. 

Remark 1.5. Notice that both and norms enhance the decay rate of the solution. The 
constraint s < 3/2 in Theorem comes from applying Lemma \A.5\ to estimate the nonlinear 
terms when doing the negative Sobolev estimates via A"^. For s > 3/2, the nonlinear estimates 
would not work. This in turn restricts p > 1 in Corollary \1.2\ by our method. However, the 
constraint p > 1 also seems necessary for the usual LP -LP' approach. In that approach, one 
should need to estimate the L^ norm of the A~^ acting on the nonlinear terms by using the 
linear decay rate p.9|) . This requires p > 1 for applying Lemma \A.5l 

Remark 1.6. Note that the L^ optimal decay rate of the higher-order spatial derivatives of the 
solution are obtained. Then the general optimal L'^ decay rates of the solution follow by the 
Sobolev interpolation (cf. Lemma \A.l\) . For instance, it follows from (|1.16p ~ (|1.17p that 



1 



QimL^ <CMt)\\}AV^Q{t)\\},<Co{l + t)--v--- and ||(n, V«>) (t)||^oo <C^{l + t) 



_3_ 
'2p , 



(1.21) 

We remark that Corollary \1.2\ not only provides an alternative approach to derive the L^-L'^ 
type of the optimal decay results but also improves the results in the usual LP-L^ approach, 
where one seems needing to assume that p is near 1 for the nonlinear system and the optimal 
decay rate of the higher-order (greater than 3) spatial derivatives of the solution is not clear. 

We will prove Theorem II. II bv the energy method that we recently developed in [9J. As there, 
we may use the linear heat equation to illustrate the main idea of this approach in advance. 
Let u{t) be the solution to the heat equation 

dtu- Au = in 
u\t=o = uo, 

Let —s<£<N. The standard energy identity of (|1.22p is 



(1.22) 



ld_ 

2di 

Integrating the above in time, we obtain 



V^u{t) 





2 






+ 




L2 





L2 



< 



2 

L2 



L2 



0. 



(1.23) 



(1.24) 



Note that the energy in (|1.23p is not bounded by the corresponding dissipation. But the crucial 
observation is that by the Sobolev interpolation the dissipation still can give some control on 
the energy: for — s < ^ < A^, by Lemma lA. 41 we interpolate to get 

i + 3 

7i„ 



V'u{t) 



< 



L2 



i+l + a 

Il2 



Combining (fT^ and (fT^ (with i 



\A-'u{t) 
—s), we obtain 



1,2 



> A" 



U0\\l2 



1 

e+s 



Vu(t) 



l+l + s 
L2 



1+ 

L2 



(1.25) 



(1.26) 
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Plugging ()1.26p into (jl.23p . we deduce that there exists a constant Co > such that 



L2 



L2 



1 

' e+s 



< 0. 



Solving this inequality directly, and by ()1.24p . we obtain the following decay result: 

-i£+s) 



L2 



< 



l+s ^ CqI 

L2 £ + s 



<Co(l + t) 



-(l+s) 



(1.27) 



(1.28) 



Hence, we conclude our decay results by the pure energy method. Although ()1.28p can be proved 
by the Fourier analysis or spectral method, the same strategy in our proof can be applied to 
nonlinear system with two essential points in the proof: (1) closing the energy estimates at each 
i-th level (referring to the order of the spatial derivatives of the solution); (2) deriving a novel 
negative Sobolev estimates for nonlinear system which requires s < 3/2 (n/2 for dimension n). 

In the rest of this paper, except that we will collect in Appendix the analytic tools which will 
be used, we will apply the energy method illustrated above to prove Theorem 11.11 However, 
we will be not able to close the energy estimates at each i-th level as the heat equation. This 
is caused by the "degenerate" dissipative structure of the NSP system when using our energy 
method. More precisely, the linear energy identity of the problem reads as: for k = 0, . . . , N , 

~ f \V''0\^ + \V''u\^ + \V''V<^fdx+ [ n\VV''uf + {fi + X)\dwV''ufdx = 0. (1.29) 
2 at J]g3 J^3 



The constraint ()1.2p implies that there exists a constant (Tq > such that 

2 



/i|VV''np + (^ + A)|divV''ti|^(ix > uq 



L2 



(1.30) 



Note that (|1.29p and (|1.30p only give the dissipative estimate for u. To rediscover the dissipa- 
tive estimate for g and V^", we will use the equations via constructing the interactive energy 
functional between u and V q to deduce 



VV^gdx + C 





2 

+ 




2 

+ 




2 

+ 


Yfc+2y^ 




L2 




L2 




L2 





< 





2 

+ 




2 




L2 




L2 



(1.31) 



This implies that to get the dissipative estimate for g and V<I> it requires us to do the energy 
estimates ()1.29p at both the fe-th and the {k + l)-th levels (referring to the order of the spatial 
derivatives of the solution). To get around this obstacle, the idea is to construct some energy 
functional £'7(t), Q<l<m-l with l<m<N {I less than iV - 1 is restricted by (fOT]) ). 

(1.32) 



£<k<in 



L2 



which has a minimum derivative count i. We will then close the energy estimates at each ^-th 
level in a weak sense by deriving the Lyapunov-type inequality (cf. (|4.4p ) for these energy 
functionals in which the corresponding dissipation (denoted by 2?™(t)) can be related to the 
energy £Y^{t) similarly as (jl.26p by the Sobolev interpolation. This can be easily established for 
the linear homogeneous problem along our analysis, however, for the nonlinear problem (|2.ip it 
is much more complicated due to the nonlinear estimates. This is the second point of this paper 
that we will extensively and carefully use the Sobolev interpolation of the Gagliardo-Nirenberg 
inequality between high-order and low-order spatial derivatives to bound the nonlinear terms 
by y^£'q (t)P^(t) that can be absorbed. When deriving the negative Sobolev estimates, we need 
to restrict that s < 3/2 in order to estimate A~'^ acting on the nonlinear terms by using the 
Hardy-Littlewood-Sobolev inequality, and also we need to separate the cases that s G (0, 1/2] 
and s G (1/2, 3/2). Once these estimates are obtained. Theorem 11.11 follows by the interpolation 
between negative and positive Sobolev norms similarly as the heat equation case. 
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To end this introduction, we will compare the NSP system (jl.ip with some related models 
with the electric force. The mostly related model is that the compressible bipolar Navier-Stokes- 
Poisson system (BNSP) of describing the dynamic of charged particles of two carrier type (e.g., 
ions and electrons). It is observed in |23[ [13] that the BNSP system can be reformulated to 
an equivalent system consisting of the NS system for the sum of densities and velocities, {p,u), 
and the NSP system for the difference, {d,w,'V^), which are coupled with each other through 
the nonlinear terms. Then for the linearized BNSP system, {p, u) decays as the NS system 
and {d, w, V<I>) decays as the NSP system. However, for the nonlinear BNSP system there is a 
new difficulty arising when estimating the nonlinear interactive terms. Precisely, there is one 
term givV^ that we can not bound it through our energy method since norm of these three 
functions are all not included in the dissipation rate. Hence, it is interesting to apply our energy 
method to the BNSP system. But in the study of the kinetic models of the Vlasov-Poisson- 
Boltzmann system, the situation is contrary. The one-species Vlasov-Poisson-Boltzmann system 
[3] will encounter the similar difficulty as the BNSP system when applying our energy method. 
But due to the special cancelation property between two species which gives the dissipative 
estimates of the norm of the electric field, in [35j we have successfully applied our energy 
method to the two-species Vlasov-Poisson-Boltzmann system to show the decay rates of the 
solution. The decay result itself is very attractive: the total density of two species of particles 
decays at the optimal algebraic rate as the Boltzmann equation, but the disparity between 
two species and the electric field decay at an exponential rate. Finally, it is also interesting to 
compare with the compressible Euler-Poisson system (EP) without the viscosity. It was shown 
in [7| that for the irrotational EP system both density and velocity decay in L°° norm at the 
rate (1 -|- for any /3 e (1, 3/2). Compared to our result (cf. ()1.2ip ) it can be comprehended 
that the dissipative effect of the viscosity enhances the decay rate of the density. 

2. Energy estimates 
Denoting g = p — 1, we rewrite (jl.ip in the perturbation form as 
dtQ + divu = —div{gu) 

dtu - nAu - (/U + A)Vdivu + V^i - V^> = -n • Vn - h{g) {fiAu + (/U + A)Vdivu) - f{g)Vg 
A$ = 

{g,u)\t=Q = {qo,uo), 

(2.1) 

where the two nonlinear functions of g are defined by 



Kg) :-- 



and f{g) := 



1. 



(2.2) 



g+1 g + l 

In this section, we will derive the a priori energy estimates for the equivalent system ()2.1 
Hence we assume a priori that for sufficiently small 5 > 0, 



'£iit) = imwrn + Mt)\\m + W'^mwm < ^- 

First of all, by ()2.3p and Sobolev's inequality, we obtain 

1/2 < g + l <2. 

Hence, we immediately have 

\hig)\, \fig)\ < C\g\ and \h^'HQ)l \f^'HQ)\ < C for any A: > 1. 
We first derive the following energy estimates which contains the dissipation estimate for u. 
Lemma 2.1. // v^f^(t) < 6, then for k = 0, . . . , N , we have 



(2.3) 
(2.4) 
(2.5) 



d 
di 



V^g\^ + \Vu\^ + I V'V^I^ dx + C 



7fc„,|2 



2 

L2 



( 




2 

+ 




2 

+ 








L2 




L2 





L2 



(2.6) 
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Proof. Applying V*^ to (j2.ip ^. ()2.iP n and multiplying the resulting identities by V^g, S/^u re- 
spectively, summing up them and then integrating over by parts, we obtain 

~ [ \V''q\^ + \V^u\^dx+ [ /i[V'=+iu|2 + (^ + A)|V'=divn|2dx- / VV''<^ ■ V'^udx 

= / V''{-0divu-u-VQ)V''g-V''{u-Vu + h{0){fiAu + {fi + X)Vdivu) + f{Q)VQ) -V^udx 
:= Ji + J2 + h + Ja + 'h- 

(2.7) 

We shall first estimate each term in the right hand side of ()2.7p . The main idea is that we 
will carefully interpolate the spatial derivatives between the higher-order derivatives and the 
lower-order derivatives to bound these nonlinear terms by the right hand side of ()2.6p . First, 
for the term Ji, employing the Leibniz formula and by Holder's inequality, we obtain 



Ji :=- [ \7^iQdivu)V''Qdx = - V / V^gV^-^divuV'' Qdx 

Jr3 



0<£<k 



0<l<k 



L2 



(2. 



L2 



To estimate the first factor in the above, we take the L°°-norm on the term with less number of 
derivatives. Hence, if ^ < [|] , together with the Sobolev interpolation of Lemma lA.H we have 



L2 



L2 



L2 



L2 



(2.9) 



L2 



Here a comes from the adjustment of the index between the energy and the dissipation, which 
is defined by 



a 1 



a 



X 1 



2,k 



2{k-£) 



+ 



:,3 



3~2j' k 
since t < —. 



(2.10) 



Hence, by the definition of the energy £q and Young's inequality, we obtain that for I < [|] 



L2 



L2 



+ 



L2 



(2.11) 



If [|] + 1 < ^ < A; (if /c < [|] + 1, then it's nothing in this case, and hereafter, etc.), we have 

l-\ 



L2 



lL2 



L2 



|V"u| 



L2 



L2 



(2.12) 



< ^1 





+ 








L2 







where a is defined by 
k + l 



a 1 
3 ~ 2 
k - 1 



1 



X 



a 



k-1 
+ 2 G 



+ 



k + l 1 



X 1 



k-1 



5,3 

2 



since i > 



k + l 



(2.13) 



2(^-1) 

In light of (|2.1ip and ()2.12p . by Cauchy's inequality, we deduce from ()2.8p that 

2 



( 




2 

+ 








L2 





L2 



(2.14) 
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Next, for the term J2, we utilize the commutator notation ()A.6p to rewrite it as 

J2-=- [ V^{u-VQ)V^Qdx = - f {u-VV^Q+[V^,u]-VQ)V^Qdx 
:= J21 + J22- 

By integrating by part, by Sobolev's inequaUty, we have 



J2i = -/ n- V ^^/^ dx = l [ divu\V''g\'^dx 
7r3 2 2 J^3 



< 



\\S/u\ 



L2 



a/ '^n 



L2 



We use the commutator estimate of Lemma I A. 31 and Sobolev's inequality to bound 



^22< llVnll^. 



2 



L2 



+ 



L2 



+ 



2 

L2 



L6 



L2 



In light of (I2J6]1 - (I217]1 . we find 

^2 < a/^:^ 





2 

+ 






L2 





We now estimate the term J3. By Holder's and Sobolev's inequalities, we obtain 



J3 = - /" V'^ (n • Vn) ■V^udx = - Y Ci I (v^u ■ VV^-'u) ■ V' 



0<f<fc 



0<£<fc 



6 



L2 



If £ < [|] , by Holder's inequality and Lemma lA.ll we have 



V^n- V 



U 



f < 



V 



L2 



1 

II V fi|| 



where a is defined by 

i_l 
3 ~ 3 



a 1 
3 ~ 2^ 



L2 



11 k + 1 



L2 



L2 



L2 



"1 



+ 



k + 1 1 



A; + 1 



, 1 since i < —. 
2' / - 2 



2(A; + 1 -£) 

If ^ > [|] + 1, by Holder's inequality and Lemma lA.ll again, we have 



V^u- V 



u 



r < 



< 



u\ 



L2 



L2 



V 



L3 

'= + 1 II V7Q., II k + 1 



L2 



||V"u 



L2 



1- 



L2 



>, A/ Cq 



L2 
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where a is defined by 

k-e + 1 I _ fa 1 
3 3 ~ 

k + 1 (\ 

=^ a -- 



3 2/ 13 2/1 



since I > 



(2.23) 



(2.24) 



In light of (ICTD and (g^l, we deduce from dXTU]) that 

Next, we estimate the term J4. We do the approximation to simphfy the presentations by 
J4:=-/ V*^'(/i(^)(/^An + (^ + A)Vdivn)) • V'^udx « - / (/i(^)V\) • V'^u dx. (2.25) 
If A: = 0, by the fact (j2.5p and Holder's and Cauchy's inequalities, we obtain 



' +iiv«iii. 



J4 « - / h{g)V^u -udx < \\q\\^2 ||V^n|| 3 ||n||^6 < a/'^o 
For /c > 1, we integrate by parts to have 

< ||V^/1(^.)V^ 



(2.26) 











L2 


IIl2 



(2.27) 



o<e<k-i 

If < £ < [|] , by (I23D, LemmaQand the estimates in (fZOjl . we obtain 





< 








L2 




< 1 




L2 





< 



L2 



L2 



(2.28) 



If [|] + 1 < ^ < /c — 1, we rewrite this factor to have 



Vh{Q)V 



L2 



{h'{0)V0) V 



L2 



0<m<£-l 



L2 



< Yl ||v™/i'(i))V^-"^V^-^+in| 



0<m<£-l 



IL2 



For m = 0, by the fact ()2.5p and the estimates in ()2.12p . we have 



h'{Q)V^QV 



u 



L2 



L2 



L2 



+ 



L2 



(2.29) 
(2.30) 



For 1 <m<^ — 1, by Lemma I A. 21 and Lemma lA.H we have 
V'"/i'(^)V^-"^V^-^+in 



<P'"'h'{e)\\^. 



L2 



.-2 „^2^„l- 



L2 

7 l|V"u||2- 



2^11"^ fc-2 
L2 



L2 
1 







L 






L2 


















L2 




U 


L2 


(2.31) 



£-1 

" k-2 



L2 





+ 








L2 
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where a is defined by 



£ -I ( £ -1 

k-i+l = ax + {k + l) X 1 



a = 3 



k-2 
k-2 



k-2 



k + 1 

^ € (1, 2] since <£<k-l. 



(2.32) 



In light of (ICTD and (I23TD . we deduce from (l2:29|) that for [|] + 1 < ^ < A: - 1, 



L2 



L2 



+ 



L2 



This together with (f2:28|) and (f2:26|) imphes that 



( 




2 

+ 








L2 





(2.33) 



(2.34) 



Finally, it remains to estimate the last term J5. If /c = 0, by the fact (j2.5p and Holder's and 
Cauchy's inequalities, we obtain 



h = - ! f{g)Vg-udx < Unlike < Js^ {mI^ + ||Vn||i.) . 

For /c > 1, we integrate by parts to have 

^5=/ v^-l(/(^>)v^.)•v'=+lndx= V C7f_i / vV(^)v'='^^)•v^+l^xc^x 

^ E IK'-^(^)^'"'42||^'^'HL.- 

o<e<k-i 

If < £ < [|] , by Lemma ED and Lemma 01 we have 



(2.35) 



(2.36) 



L2 



VV(£>) 



L2 



< \N>\\p 



V^g 



L2 



ll^'ll 



1-1 
L2 



y '-'0 



L2 



(2.37) 



L2 



where a is the same one defined by (I2TI1 . If [|] + 1 < ^ < A; - 1, we rewrite this factor to have 



v'f{0)^'-'0 , = ||v^-i (/'(^)v^) v^-^^, 



L2 



o<m<e-i 



L2 



0<m<^-l 



For m = 0, by the fact (|2.5p and Lemma lA.li we have 
f'i0)V'gV''-'g 



< 








L2 ~ 




L2 





V^g 



l|V"£>|||. 



1- 



L2 



V^g 



(2.38) 



(2.39) 



L2 



where a is defined by 



DECAY OF NSP 

For 1 < m < i — 1, hy Lemma I A. 21 and Lemma I A. 11 we have 
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< ||V"'/'(e)|L 



L2 



^ II 9 i|2-- 



,1- 



L2 



V'-'^g 



V'-'g 



L2 



L2 " ' 



V'^g 



lL2 

''-^ 'IV^nll'-'"^ 



V'^g 



L2 



L2 



L2 



V'^g 



L2 



~ ^0 



L2 



where a is defined by 

k — i = a X 



1 , ^-1 



/c-2 
k-2 



e (0, 11 since < ^ < A; - 1. 



L2 



>. A/ Cq 



L2 



In light of and dXIT]) . we deduce from (CT]) that for [|] + 1 < £ < A: - 1, 

V'fig)V'^-'g 
This together with ()2.37p and (|2.35p imphes that 

Now we turn to estimate the left hand side of (|2.7p . For the second term, we have 
/ fi\V''^^u\'^ + {fi + X)\V''divu\'^dx>ao 



( 




2 

+ 








L2 





2 

L2 



(2.41) 



(2.42) 



(2.43) 



(2.44) 



(2.45) 



While for the last term, by the continuity equation ()2.ip ^ and the Poisson equation ()2.ip o and 
the integration by parts, we get 



VV''^ ■V'^udx = / V'^^V'^divudx 



1 d 
2di 



[ V''<^V''dtg + V''<^>V^dw{gu)dx 
[ V'=$V'=5i A$ - V^'(^^) • V'V^* dx 
[ |V'=V$pdx+ / V^gu) ■V'^V^dx. 



(2.46) 



Notice carefully that we can not estimate the last term in (|2.46p directly. For instance, we may 
fail to bound f^^ gV^u-V'^'V^ dx by the right hand side of (|2.6p . To overcome this obstacle, the 
key point is to make full use of again the Poisson equation ()2.1P 3 to rewrite g = A<I>. This idea 
was also used in [31]. Indeed, using ()2.ip o and the integration by parts, by Holder's inequality 
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and Lemma lA.ll we obtain 

Jr3 J]r3 

= / V^(V^> • Vu) ■ V^'V$ + V'=(V$ • u) ■ VV^V$ dx 



0<£<k 



(2.47) 



o<e<k 



o<e<k 

















L3 









+ 












L2 




LB 




LB 







L2 



L2 



If ^ < [|], by Lemma EH we have 









^ U 











L2 



< ||V°V<I>| 



L2 



t 

k + 1 
L2 



I fe+i 
Il2 



1- 



fc+i 



L2 



>, A/ Cq 



( 




+ 










L2 







(2.48) 



where a is the same one defined by ([2:21]) . Now if £ > [|] + 1, by Lemma [O 

again, we have 
(2.49) 

















L2 



f+1 



^2 I' IIL2 



fc+i 



L2 



^ V ^0^ ( 





+ 








L2 







where a is defined by 



/c — £ + l = ax 



+ 



k + l 



+ {k + l) 



+ 



fc + 1 
a = r S 



1 \ , /c+l 

-, 1 I smce £ > . 



(2.50) 



2£ + l 

In hght of the estimates ([2:i8]) ~ (f2:i9]) . we deduce from ([Tie]) that 

- 2dt 



( 




2 

+ 








L2 





L2 



(2.51) 



Plugging the estimates for Ji ~ J5, i.e., QUI) . (|2J8]) . ([2231), (1231) and ([233]), and the 
estimates (l23T]) and ([TIS]) into (H^D, we get ([Ml). □ 

The following lemma provides the dissipation estimate for q and V$. 

Lemma 2.2. // y^oW ^ ^^^^ A; = 0, . . . , — 1, we have 



d 
lit 



V^u-VV^gdx + C 





2 

+ 




2 

+ 




2 

+ 


Yfc+2y^ 




L2 




L2 




L2 





< 





2 

+ 




2 




L2 




L2 



(2.52) 



Proof. Applying V'^ to (|2.1|) a and then taking the inner product with W^g, we obtain 
/ \VV^Q\'^dx- [ V^V<I> ■VV'^Qdx 



< 



W''dtu-VV''Qdx + C 











L2 





L2 



(2.53) 



+ 



V'' (n • Vn + h{Q){nAu + {p + A)Vdivu) + f{Q)Vg) 







L2 





L2 
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The delicate first term in the right hand side of (j2.53p involves the time derivative, and the 
key idea is to integrate by parts in the t-variable and use the continuity equation. Thus by 
(j2.ip ^ and integrating by parts for both the t- and x- variables, we may estimate 



-4- I y'"u-VV^Qdx- 
at ./iii3 



V^div-u • V^Qt dx 



(2.54) 



V^'divu 



L2 



+ / Vdivn • VdiwiQu) dx. 



By Holder's inequality, we have 

/ V^div-u • vMiv(^n) < ^ ||v^^V^+^- 
•^^^ 0<£<fc+l 



L2 



L2 



(2.55) 



If < £ < [^] , by Lemma [AH we have 



V'qV 



u 



L2 



V 



fc+i- 



I 

k + 1 



L2 

e 

k+i 
L2 



I k + 1 

Il2 



fe+1 



L2 



(2.56) 





+ 








L2 







where a is defined by 



a 1 
3 ~ 2 



1 



k + 1 



+ 



a 



Sjk + l) 

2{k + i-e) 



< 3 since £ < 



k + 1 1 
~~3 2 
k + 2 



k + 1 



(2.57) 



While for £ > [^y^] + 1 (then k + 1 — £ < [^y^] ), we can then interchange the roles of g and u 
to deduce that (|2.56p holds also for this case. Thus, in view of (|2.54p - (|2.56p . we obtain 



V-Ui • W'gdx < -— / ■ VV^gdx + C 

dt 



L2 



L2 



(2.58) 



Next, note that it has been already proved along the proof of Lemma |2. II that 
V'' {u-Vu + h{g){fiAu +{fi + A)Vdivu) + f{g)S/g) 



L2 





+ 








L2 







(2.59) 



We now use the integration by parts and the Poisson equation (j2.iP o to have 
V^V<I> ■W'^gdx = [ \/^A^V''gdx= [ \V''g\'^dx. 



On the other hand, it follows from the Poisson equation that 



L2 



L2 



L2 



and 



L2 



V'^+^g 



L2 



(2.60) 



(2.61) 



Consequently, by (|2.58p -- (|2.6ip . together with Cauchy's inequality, since \/£q < 5 is small, 
we then deduce (p32|) from ([233]) . □ 



14 yanjin wang 

3. Negative Sobolev estimates 

In this section, we will derive the evolution of the negative Sobolev norms of the solution to 
(j2.ip . In order to estimate the nonlinear terms, we need to restrict ourselves to that s € (0, 3/2). 
We will establish the following lemma. 

Lemma 3.1. // y^fp (t) < 6, then for s € (0, 1/2], we have 



dt Jr3 " 

^ (II^IIh^ + ||Vn||^i) + \\^''u\\l2 + llA^'V^ll^a) ; 

and for s & (1/2,3/2), we have 

dt Ju3 " 



(3.1) 



(3.2) 



Proof. Applying A to ()2.1P ]^. ()2.ip n and multiplying the resulting identities by A A '^u 
respectively, summing up them and then integrating over by parts, we obtain 

~[ \A-'Q\'^ + \A-'u\'^dx+ [ /u|VA-"n|2 + (^ + A)|divA-"u|2dx- / A^'V^ ■ A^'udx 

= / A""* (-gidivu - u ■ Vg) A'^'g - A'" (n • Vn + h{Q){i^Au + (/x + A)Vdivu) + f{g)Vg) ■ A^^udx 
:= Wi + W2 + W3 + Wi + W5. 

(3.3) 



In order to estimate the nonlinear terms in the right hand side of ()3.3p . we shall use the 
estimate ()A.13p of Riesz potential in Lemma lA.51 This forces us to require that s G (0, 3/2). If 
s G (0, 1/2], then l/2 + s/3 < 1 and 3/s > 6. Then by Lemma lA.51 and Lemma lA. 11 together 
with Holder's and Young's inequalities, we have 



Wi = - [ A-'{Qdivu)A-'gdx < l|A-"(gdivw)L2 ||A-"g| 

7R3 



L2 



< ||^divn||^_^_i_ ||A "^11^2 < \\g\\L3/s ||Vu||^2 ||A "^||^2 

< ||V^||/2 ||V £']|^2 IIV-ull^ajlA q\\^2 



(3.4) 

\\^''e\\L2 

Similarly, we can bound the remaining terms by 

W2 = - j^^ A-\u ■ Vg)A~'gdx < (||Vn||^i + WVgWl^) \\A-'q\\l2 , (3.5) 

W3 = -J A^'{u-Vu)-A-%dx<(^\\Vufj^^ + \\\/^u\\l2)\\A"'u\\^2, (3-6) 

W4 = - [ A-' {h{g){nAu+ + X)Vdivu))A-'udx (3.7) 

< (l|Vg||^i + ||v2^x||^2) ||A-^^x||^2, (3.8) 

W5 = -J^^ A-« {f{g)Vg) ■ A'^udx < [\\S/g\\l^ + WV^gWl^) \\^''4l^ ■ (3-9) 

Now if s € (1/2,3/2), we shall estimate the right hand side of ()3.3p in a different way. Since 
s G (1/2,3/2), we have that 1/2 + s/3 < 1 and 2 < 3/s < 6. Then by Lemma I A. 51 and Lemma 
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lA.ll we obtain 



■ / A ''{gdi-vu)A ''gdx<\\A ''(^divu)|| 2 ||A ""9^2 



|Vu||^2 IIA""*^ 
Similarly, we can bound the remaining terms by 

W2 



A '^11^2 < I|Vu||^2 ||a 'q 



IL2 



(3.10) 



^ II iiS-1/2 IIV7 i|3/2-s 
< \\q\\l2 



L2 



/ A~'iu ■ VQ)A-'Qdx < \\u\\'j2^^^ ||Vn||5^/^"' ||V^||^2 ||A-"^| 



L2 



A-" {u ■ Vu) ■ A-'udx < Ml-,^^^ ||Vu||5^( 



2-s 



Vn||^2 ||A ""ul 



L2 



Wi = - i A"^(/i(^)(^An+ (/i + A)Vdivu))A"^udx 



^ II ll«-l/2 



|3/2-s 
Il2 



\A-'u\ 



(3.11) 
(3.12) 

Il2 11^'^ "'1Il2 ' (3.13) 

' |Vi?||^2||A-^^||^2- (3.14) 

Finally, we turn to the left hand side of (|3.3p . For the second term, we have 

/ n\VA-'u\^ + {li + \)\d{^iA-'u\^dx- [ A-'V<i>- A-'udx>ao\\VA-'u\\%. (3.15) 
Jm.3 J^3 

While for the Poisson term, by the continuity equation ()2.ip ^ and the Poisson equation ()2.ip n 
and the integration by parts, we get 

A-'^V^ • A-'udx 



[ A-' {f{Q)VQ) ■ A-'udx < WeWl-'^' \\Vq\ 



/ 



A *$A ^divudx 



[ -A~'^A~'dtQ - A-"^>A-Miv(^n) dx 
Jrs 

/ -A~'<^A-'dtA<i> + A-'V<i>-A-'{Qu)dx 

JR3 

— ^ / \A-'V^\'^dx+ I A-'V^ ■ A-'{gu)dx. 
21 dt Jks 



(3.16) 



If s G (0, 1/2], we use Lemma lA.51 and Lemma lA. II to obtain 

IIa-s/ \II ^11 ^ x-^ II 1/2— s ||v-^2 l|l/2+« 

||A (£''")||i2 < bllL2 |ln||^3/. < 11^11^2 llVnll^'a ||V n|i 



and if s € (1/2,3/2), we have 



|A '{qu)\\^2 < I|£'IIl2 hh^/s < I|£'IIl2 lkllz,2 



s-1/2 



Vn 



IL2 



|3/2-s 

Il2 ' 



Consequently, in hght of ([HI^ - ffXTS]) . we deduce ([321) from l^^i). 



(3.17) 

(3.18) 

□ 



4. Proof of Theorem HJJ 



In this section, we shall combine all the energy estimates that we have derived in the previous 

two sections and the Sobolev interpolation to prove Theorem 11.11 

We first close the energy estimates at each £-th level in our weak sense to prove ()1.12p . Let 
> 3 and < I < m — 1 with 1 < m < A. Summing up the estimates (|2.6p of Lemma 12.11 for 

from k = i to m, since \/£q < 5 is small, we obtain 



d_ 
di 



e<k<m 



L2 



+ 



L2 



+ 



e+l<k<m+l 



2 

L2 



,£<k<m 



i+l<k<m+l 



2 

L2 



(4.1) 
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Summing up the estimates (j2.52p of Lemma 12.21 for from k = i to m — 1, we have 



I E X,v'„.vvV. + cJi: ||v' 

£<fc<m-l'''* \i<k<m 



L2 



+ 



E 



£+l<k<'m+l 



v^v$ 



2 

L2 



(4.2) 



e+l<k<m+l 



L2 



Multiplying (|4.2|) by 2C2(5/C3, adding the resulting inequality with (|4.ip . since (5 > is small, 
we deduce that there exists a constant C5 > such that for < £ < m — 1, 



E 

l<k<m 





2 

+ 




2 

+ 






L2 




L2 





l+l<k<m+l 



i<k<m 



2 

L2 
2 

L2 



e<k<m-l 



(4.3) 



^+l<A;<m+l 



2 

L2 



< 0. 



We define (t) to be C^^ times the expression under the time derivative in ()4.3p . Observe 
that since 5 is smah, ^7^(i) is equivalent to ||V^^(Ot„-£ + || V^n(t)||;^„_, + || V^V$(t)||^„_,. 
Then we may write ()4.3p as that for < £ < m — 1, 



+ 




2 

+ 




2 

+ 








^m — l 









< 0. 



(4.4) 



Now taking £ = and m = 3 in ()4.4p and then integrating directly in time, we get 

Mml + \Ht)\\l^ + WvmWl'^ < slit) < £1(0) < WqoWI + WmWl. + ||v<i>o||^3 . (4.5) 

By a standard continuity argument, this closes the a priori estimates (|2.3p if at the initial time 



we assume that H^^olls + ||'Uo|Ih3 + 1 1 ^'^o 1 1^3 — "^o is sufficiently small. This in turn allows us to 
take i = and m = in (|4.4p . and then integrate it directly in time to obtain (jl.l2p . 

Next, we turn to prove ()1.13p - (|1.15p . However, we are not able to prove them for all s € 
[0, 3/2) at this moment. We shall first prove them for s € [0, 1/2]. 

Proof of (frT3]l - (frT5D forse [0,1/2]. Definef_,(t) := ||A-^^(t)||^2 + ||A-"u(t)||^2 + ||A""V$(t)|| 
Then, integrating in time (|3.ip . by the bound (|1.12p . we obtain that for s S (0, 1/2], 

£-s(.t) < £-s{0) + C [ (WqWI^ + ||Vu|]^i) ^£-s{r)dT 

^° ^ ^ (4.6) 



< Co ( 1 + sup V£-siT) ] . 

0<T<t 



This implies (fLTHIl for s e [0, 1/2], that is. 



\\A~'Q{t)\\l, + \\A-'u{t)\\l, + \\A-'V^t)\\l, < Co for s G [0,1/2]. 
If ^ = 1, . . . , — 1, we may use Lemma lA. 41 to have 



By this fact and (j4.7p . we may find 



L2 



> c||A-7||^/+^ 



L2 



+ 



L2 



L2 



+ 



L2 



L2 



(4.7) 



(4.8) 



(4.9) 
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This together with (jl.l2p imphes in particular that for i = 1, . . . , N — 1, 

2 





2 




2 








+ 




HN-e + 









> Cn 





2 




2 










+ 















jjN-e-i 
2 



(4.10) 



Thus, we deduce from (|4.4|) with m = N the following time differential inequality 

,N-1. 



+ Co(^f)'+^ <0 for^ = l, 



Solving this inequality directly gives 

£^{t) < Co(l + for £ = 1, . . . , 4iV - 1. 

This implies that for s € [0, 1/2], 



(4.11) 



(4.12) 



+ 



HN- 



+ 



vV«>(t) 



ffN+i-e 



On the other hand, since g = divV<I>, we have 



< Coil + t)-(^+^) for £ = 1, . . . , TV - 1. 

(4.13) 



V^£.(i) 



L2 



< 



V^+iV$(t) <Co(l + t)"(^+^+") for 



0,...,N-2. (4.14) 

Hence, by (H33l) . KIM , (g^]) and the interpolation, we get (fna - (frT5]l for s € [0, 1/2]. □ 

Now we can present the 

Proof of (frT3]) - (frT5|) for s G (1/2,3/2). Notice that the arguments for the case s £ [0, 1/2] 
can not be applied to this case. However, observing that we have go,uo,V^o G H~^/^ since 
H^'^ n C for any s' € [0, s], we then deduce from what we have proved for ()1.13p - (|1.15p 
with s = 1/2 that the following decay result holds: 



and 



+ 



+ 



vV$(t) 



ffN+i-e 



^^Qit) < Coil + t)-^^+^/^'^ for £ = 0,..., TV -2. 



< Co(l + t)"(^+^/^) for £ = 0,...,iV-l. 

(4.15) 

(4.16) 



Hence, by ()i35]) ^ (fmD . we deduce from ([32]) that for s G (1/2,3/2), 

£.s{t) < £-s{0) + C f \\{g,u)r^,^'^ {\\g\\H2 + \\Vu\\H.f^-' V£ZF)dT 
Jo 

<Co + Co [\l + T)-^"^-''^Ur sup ^/sIJt) 

Jo 0<T<t 

<Co(l + sup VsZir)) ■ 

V 0<T<t / 

This implies (fTTHD for s G (1/2,3/2), that is, 

\\A-'g{t)\\l, + \\A-'u{t)\\l, + \\A-'Vm\\l2 < Co for s G (1/2,3/2). 



(4.17) 



(4.18) 



Now that we have proved (|4.18p . we may repeat the arguments leading to (|1.14p - (|1.15p for 
s G [0, 1/2] to prove that they hold also for s G (1/2, 3/2). □ 
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Appendix A. Analytic tools 

A.l. Sobolev type inequalities. We will extensively use the Sobolev interpolation of the 
Gagliardo-Nirenberg inequality. 



Lemma A.l. Let < m, a < ^, then we have 

l|V"/||^. < ||v-/||i7^ 

where a satisfies 



a _ 1 
3 



m 1 



1 



(A.l) 

(A.2) 

□ 



Proof. This is a special case of pHl pp. 125, THEOREM]. 

Next, to estimate the L°° norm of the spatial derivatives of h and / defined by (|2.5p . we shall 
record the following estimate: 

Lemma A.2. Assume that \\Q\\i^aD < 1. Let g{Q) he a smooth function of g with bounded 
derivatives of any order, then for any integer m > 1 we have 



l|V'"(5(f?))ILoo < ||V>|L^. 

Proof. Notice that for m > 1, 

V™(5(^)) = a sum of products g'^''-'^"{g)V^' q ■ ■ ■ V^"^, 

where the functions (^'''^'•••''^" (£») are some derivatives of g{g) and 1 < 7i < m, i = 1, 
7i + • • • + 7n = ""T-. We then use the Sobolev interpolation of Lemma lA.ll to bound 

l|V'"(<7(£^))!Loo<||V^VllLoo---||V^"^>||i^ 

< lUII^^'^i/"^ \\\7'^n\n'^/"^ lUII^^"^"/'" IIW^^II'''"/'^ < 

rs^ 1 1 t: 1 1 ]^'^ 

Hence, we conclude our lemma since ||^>||^oo < 1- 
We recall the following commutator estimate: 
Lemma A. 3. Let m > 1 be an integer and define the commutator 

[V^J]g = V'^{fg)-fV'^g. 

Then we have 

II[V'",/]5|Il. < l|V/||i.. ||V"-^<7|L.. + l|V"^/||^.3 



|V"^^|| 



where p,p2,P3 G (l,+oo) and 



1 _ 1 1 _ 1 1 

P Pi P2 P3 Pi 



(A.3) 

(A.l) 
, n with 

(A.5) 

□ 



(A.6) 
(A.7) 
(A.8) 



Proof. For p = p2 = = 2, it can be proved by using Lemma lA.li For the general cases, one 
may refer to [14, Lemma 3.1] □ 



A.2. Negative Sobolev norms. We define the operator A^,s € M by 



A7(x) 



(A.9) 



where / is the Fourier transform of /. We define the homogeneous Sobolev space of all / 



for which 



is finite, where 



:= ||A^/||^. = \i\^f 



L2 



(A.IO) 



We will use the non-positive index s. For convenience, we will change the index to be "— s" 
with s > 0. We will employ the following special Sobolev interpolation: 



Lemma A. 4. Let s > and ^ > 0, then we have 



< 



L2 



1- 



L2 



\f\\%-s , where 



+ 1+S 



(A.ll) 
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Proof. By the Parseval theorem, the definition of (jA.lOp and Holder's inequaUty, we have 



L2 



< 



L2 



1-6 



L2 



L2 



1-6 



L2 



(A.12) 

□ 



If s G (0,3), A defined by (|A.9p is the Riesz potential. The Hardy-Littlewood-Sobolev 
theorem implies the following type inequality for the Riesz potential: 



Lemma A. 5. Let 0<s<3, 1 < p < q < oo, 1/q + s/3 = 1/p, then 



\A-'f\\r.<\\fh.- 



Proof. See [30', pp. 119, Theorem 1]. 



(A.13) 

□ 
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